A magnetic monopole in topological insulator: exact solution and Witten effect 
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The Witten effect tells that a unit magnetic monopole can bind a half elementary charge in an 
axion media. We present an exact solution of a magnetic monopole in a topological insulator that 
was proposed to be an axion media recently. It is found that a magnetic monopole can induce one 
zero energy state bound to it and one surface state of zero energy. The two states are quite robust, 
but the degeneracy can be removed by external fields. For a finite size system, the interference 
of two states may lift the degeneracy, and the resulting states have one half near the origin and 
another half around the surface, which realizes the Witten effect. However, the energy difference 
decays exponentially with the size of the system. The exact solution does not fully support the 
realization of the Witten effect in a topological insulator. 
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Topological insulators are electronic materials that be- 
have like insulators or semiconductors in the bulk, but 
are surrounded by a topologically protected conducting 
layer near the surface of the materials^— . The mate- 
rials have been studied extensively. One of the pre- 
dicted features in the materials is the so-called "axion 
electrodynamics"—, as an electromagnetic response to ex- 
ternal fields. The idea of the axion was first introduced to 
address the strong charge-parity problem in the physics 
of strong interaction^.. It becomes a possible candidate 
for the dark matter in the universe, and however, has not 
been yet confirmed experimentally so far. Historically, it 
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was known that an additional term Ot^-tII ■ E can be 
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introduced into the Maxwell Lagrangian for electric and 
magnetic fields, which is time reversal invariant when 
9 = tt. The additional term revises both the Gauss' law 
and Ampere's law in the Maxwell's equations by adding 
extra terms^, 



V • D = Pe - 



vT/iQC 



-V0-B, 



(1) 



V x H = dtD + j 



{W9 x E + d t 0B) (2) 



where D = e E + P and H 



ture constant a = 



fJ-o 



M. The fine struc- 



2eohc ' 



One of the fundamental prop- 
erties in the revised Maxwell's or axion equations is the 
Witten effect, which states that a magnetic monopole of 
unit strength ejv/ — 4>o — m an axion media must bind 
an electric charge, — (n+ j-)e, where e(> 0) the elemen- 
tary charge and n is an integer. Consider a point-like 
magnetic monopole situated at the origin. The magnetic 
field is given by V • B = cf>o5(r). We suppose that 9 = 
initially and then increases adiabatically to 9 — tt. 9 is 
uniform in the space and there is no current in the media. 
It follows from the axion equations that 



Sp e = p e {9 = tt) - p e 
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As a magnetic monopole cannot induce a half elemen- 
tary charge in a conventional media of 9 = 0, the charge 



bound to the monopole should be — e/2 except for an in- 
teger for a time reversal invariant topological insulator of 
9 = TT 7 . Charge fractionalization in condensed matters 
was extensively discussed in 1980s&£. The quasi-particles 
in the fractional quantum Hall effect carry fractional 
charge 10,11 . Whether or not a half elementary charge 
bound by a magnetic monopole could exist in a topolog- 
ical insulator becomes a subtle issue to test the validity 
of the axion theory for topological insulators. Rosenberg 
and Franz studied the Witten effect in a crystalline topo- 
logical insulator numericallyJ^, and intended to use it as 
a criterion to justify whether the system is topologically 
trivial or non-trivial^. 

In this paper, we present an exact solution of a mag- 
netic monopole located in the center of a sphere of a 
topological insulator, which is described by a modified 
Dirac-like equation. It is found that there exist two solu- 
tions of zero energy: one is located in the vicinity of the 
magnetic monopole and the other around the surface, 
which is characteristic of non-trivial topological insula- 
tor. At half filling, only one electron occupies the two 
degenerate states due to the particle-hole symmetry in 
the system. The double degeneracy of the zero energy 
states does not favor or disfavor the Witten effect be- 
cause the bound charge near the monopole can be from 
zero to — e, although it does not exclude a half elementary 
charge as the Witten effect requires. The degeneracy of 
the two states can be removed due to the finite size effect 
for a small sphere. In the case the occupied state have 
one half in the vicinity of the monopole, and another half 
is distributed around the surface, in which the Witten ef- 
fect can be realized. In general, the external fields such 
as the Zeeman field and electric potential can remove the 
degeneracy of the two states, but could not mix them to 
realize the Witten effect. 

The model Hamiltonian for a magnetic monopole in 
the modified Dirac equation is given by^ 
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where 2mv 2 is the energy gap between the conduction 
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band and valence band, v is the effective velocity and B 
is a parameter of dimension of inverse mass. a xy z are 
the Pauli matrices. The canonical momentum operator 
II = -ihV + eA and V x A = ||</>or/r 3 . q is half of 
an integer due to the quantization of magnetic charge 
and q = § for a unit monopole^. It is known that the 
vector potential A cannot be written as a single expres- 
sion in the whole space, and has to be defined as two 
functions in two overlapping regions to keep singular 17 . 
In the absence of the magnetic monopole, the topologi- 
cal properties of this equation have been studied exten- 
sively. It is topologically non-trivial for mB > and 
trivial for mB < 0. A topological quantum phase tran- 
sition occurs at mB = 14 i 15 . In the presence of the 
magnetic monopole, the orbital angular momentum is 
modified to L = r x II — qh-, which satisfies the alge- 
bra [L Q ,L^] = ihe a p 1 L 1 . The eigenfunctions of L 2 and 
h z are denoted by Y q ,i,i^ with the eigenvalues 1(1 + l)h 2 
and l z h. Since the two terms in L are orthogonal to each 
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The total 



§ a. The total 



other, L 2 = |r x n| 2 + q 2 h 2 and 1(1 
angular momentum is defined as J = L 
angular momentum J 2 and its z-component J 2 commute 
with the Hamiltonian, and are good quantum numbers. 
Thus we can diagonalize simultaneously H, J 2 and 3 Z . 
According to the definition, J is the sum of two angu- 



Thus the eigenvalues of J 2 
l 



lar momenta, L and er. 

be j = I + ^and j = I — ± respectively, and the corre- 
sponding eigenstates are the linear combination of Y q jj z 
and the eigenvector for a z . For a minimal l(= \q\) and 
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the eigenstates are 



Vj,. 
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kl- 



2M+i 2 Y q,\q\,iz + k 



(5) 



Due to the particle-hole symmetry in the model Hamil- 
tonian, the eigenvalues of ±E appear in pairs. For a 
half-filled system, we are interested in the energy lev- 
els near the Fermi surface at E = if it exists. Thus 
we focus on the case of I — \q\ and j = \q\ — | while 
other cases of j and I may have the eigenstates be- 
low or above than the Fermi level. To solve the prob- 
lem, using the relations^, (a ■ r) T]jj z — rsgn(q)r)jj s and 
a ■ I !./'('• :■//,. , = -isgn(q)(d r + r -1 )/(r)7frj, , where f(r) 
are arbitrary function of r, we construct a simultaneous 
eigenstate for iJ, J 2 and 3 Z : 
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The equation for the radial part of the wave function is 
reduced into 



1 + sgn(mB) 
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Here p = kr (k 2 = \mv 2 /BH 2 \), £ = sgn(qmv) / yJ\mB\ 
and A = E/mv 2 . 

We first consider a large radius limit of the sphere, 
kR 3> 1. In this case it is reduced to a one-dimensional 



modified Dirac equation by ignoring the term of \ at the 
end of r — R: the sign of mB determines whether there 
exits a solution of zero energy solution near r = Rf^. For 
A = 0, there is a general solution 



Vj,m = Xs®Vj,jJs(p) 
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( x T s =73(1,") and ( a = ±1)). For 



f s (p) = C ie -«>/ 2 ^-J a [fip) + C 2 e-<fl 2 ^-K a (pp) . 

Vp 
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where J a (x) and K a (x) (a — \J\q\ + \ and f3 — 

y/T— C 2 /4) are the first and second Bessel functions, 
and Ci (i = 1,2,3,4) is the normalization constant. For 
C 2 > 4 or an imaginary f3 the solution is still valid. We 
can also use the modified Bessel functions to replace the 
Bessel functions. Consider the boundary condition at 
p = and p = kR = pp<. We have two solutions of zero 
energy. For s — 1 by taking > without losing gen- 
erality, one has a solution which is convergent at p = 0, 
but decays exponentially with a larger p 



f + (p)^ Cl e^/ 2 ^-J a ((3p). 
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(10) 



Thus the wave function denoted by ipj~ m (p) — X+ ® 
ri m f+(p) is mainly located near the origin or the mag- 
netic monopole. On the other hand, both J a (j3p) and 
K a (j3p) become divergent for a large p, but convergent 
for a small p. For s = — 1 one has the other solution 
which vanishes at pr = kR 



f-(p) = Ci 
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Kl*/ 2 / J a (fip) K a (Hp) 
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This solution denoted by ip7 m (p) = X- ® Vmf-(p) is 
distributed near the surface of r = R and decays expo- 
nentially in pr ; — p or R — r. This is a surface state of 
zero energy, and is one of the characteristics of topolog- 
ical insulators. For C, 2 — 4 and (3 — 0, we also have two 
solutions of zero energy: one is near the origin, 



f+(p) = C 3 e 



-\Z\p/2 cc-l/2 



with s = 1 and the other is around the surface, 



f-(p) = C, 
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= 



with s = — 1. 

Except for the two solutions of zero energy, there also 
exist other states even for I = \q\ and j = \q\ — \. We solve 
Eq. ([7]) numerically, and present the results in Fig. 1. 
We note that the solutions of zero energy appear only for 
the case of mB > 0. There is no solution of zero energy 
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FIG. 1: Energy spectrum of the bound states of I = \q\ = | 
and j — I — j as a function of the dimensionless parameter 
mB. The energy unit is |m| v 2 . In general, the states of the 
zero energy is 4 |g|-fold degeneracy. 
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FIG. 2: The energy splitting of the two bound states as a 
function of the spherical radius pa = kR. The energy unit is 
\m\ v 2 and k 2 = \mv 2 / Bh 2 \. 



for mB < as the system is topologically trivial. We find 
that all other states are not the surface states, but the 
bound states like those with different total quantum num- 
bers in a hydrogen atom. As the results are independent 
of j z , all these states are 2(2j + I) = 4|g|-fold degener- 
ate. For q = ±i, the states are singlet with total angular 
momentum j = 0. Due to the topology of the band struc- 
tures in topological insulator the existence of the surface 
state is robust against the continuous deformation of the 
surface shape, impurities or interaction while the zero en- 
ergy state around the magnetic monopole is located very 
close to the origin, although the energy eigenvalues may 
shift away from the zero energy. 

In the present case of half filling, the Fermi level is 
located at E = due to the particle-hole symmetry in 
Eq. (fj| , below which all states or levels are fully filled. Al- 



though we do not obtain general solutions of other states 
of different j and j z , some of which are the surface states 
between the energy gap, the double degenerated states 
at E = are closely related to possible realization of the 
Witten effect. For q = i, there is only one electron to 
occupy these two well-separated states. When mB < 0, 
there is no solutions of zero energy near the center and 
around the surface. Although it is possible that the mag- 
netic monopole can bind a lot of electrons with the energy 
below zero, the electron charges accumulated around it 
must be an integer multiple of the elementary charge — e. 
This is consistent with the fact that the system is topo- 
logically trivial with 8 = 0. When mB > 0, the system 
is topologically non-trivial, the appearance of the sur- 
face state of E = is one of the characteristics. If the 
topological insulator is really an axion media, 9 should 
be 7T. Thus the sign change of mB should accompany the 
change of 9 from to 7r. Therefore if a topological insu- 
lator is really an axion media, a half elementary charge 
must be bound to the magnetic monopole as the Witten 
effect requires. Since the doubly degenerate states are 
well separated in space, our exact solutions do not favor 
or disfavor the picture of the Witten effect as the the 
charge bound to the monopole can be from to — e. 

Now we come to consider a sphere of a finite radius R. 
Denote the state near the origin by i^t m and the surface 
state by 4>J m ■ If the radius is large enough such that the 
two wavefunctions have no overlap in space, ^ m and 
tf>7 m are two exact solutions of zero energy. However, 
when the radius R is finite and the two wavefunctions 
overlap in space, the two wavefunctions interference with 
each other, and the degeneracy can be lifted, which is 
dubbed the finite size effect^. As an degenerated per- 
turbation approach, we still use the two functions as the 
basis. As Xs are the eigenstates of a y , we have the rela- 
tions x\ a zXs = x\ a xXs = 0. Therefore the expectation 
values (ip^ m \H\^ m ) = 0. However, xt^zX-s = 2 and 
xl^xX-s = -2is, then A = (ip^ m \H\^J m ) ^ 0. As a 
consequence of the first-order degenerate perturbation, 
the energy eigenvalues become ±|A| and the two states 
become 

^ = 71 V^ m ± W\ ^ m ) ' (14) 

The values of the ±| A | are evaluated numerically and 
plotted in Fig. 2. There are two different cases. For 
q 2 < 4, the gap increases monotonically with decreas- 
ing R while for <j 2 > 4 the gap oscillates, but the am- 
plitude increases with decreasing R. The states tp± are 
separated into two halves as an d ^Jm are orthog- 
onal and well separated in space. The weights of these 
two parts are exactly equal to 1/2. When the state with 
lower energy is occupied, the electron will be split into 
two parts: one half is near the origin and the other is 
around the surface as shown in Fig. 3. Thus the frac- 
tional charge is indeed realized in this case. To estimate 
the size of the wave package in real materials, we adopt 
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the fitted parameters from ARPES data for Bi2Se3 20 : 
mv 2 = 0.126eV, Bh 2 = 21.8eVA 2 , and hv = 2.94eVA. 
As the wave function decays exponentially combining 
a power law decay from the origin, exp[— r/£] or from 
the surface, exp[— (R — the characteristic length is 

£ ~ 15 A. The finite size effect becomes obvious when R 
is comparable with several 2£ ~ 30A, which is consistent 
with the fact that gap opening in the Bi-zSes thin film 
with thickness several quintuple layer o 15 ' 20 . Therefore 
the characteristic length for a finite size effect is quite 
small. Rosenberg and Franz claimed that a half elec- 
tron charge is really bound to a magnetic monopole in 
a crystalline topological insulator 12 . In their numerical 
calculation, the lattice size is of 20 3 for a cubic lattice. 
The finite size effect is obviously seen from the distribu- 
tion of excess charge: the size for the bound state around 
the monopole is about 4-5 lattice spaces and that for the 
surface state is about 6-7 lattice spaces. Thus what ob- 
served in their calculation is actually a finite size effect. 
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FIG. 3: The radial probability density Y{p) = p 2 (|/| 2 + \g\ 2 ) 
with a radius pn = kR = 60. The energy unit is |m|n 2 and 
k 2 = \mv 2 /Bh 2 \. 



Another question is whether charge fractionalization 
can be realized by the Zeeman field in the topological 
insulator or due to the proximity effect of a ferromagnetic 
layer covering around the surface. Consider a radial field 
B = B r ~. The Zeeman energy term is 



SH = fi e ffB r 







(15) 



which breaks the time reversal symmetry. The energy 
shifts of the two states of zero energy are AE± — 
2sgn(q)(j, e ff Jq R dpB r \f±(p)\ 2 . The degeneracy will be 
removed as the two integrals are usually not equal. For 
example, if the field appears only near the surface due to 
the proximity effect of a ferromagnetic layer, it is obvi- 
ous that AE + — and AE- ^ 0. However, as the off- 
diagonal integral is zero, (i>j m \ SH \ i>7 m ) = 0, the two 
states will never be mixed as those due to the finite size 
effect. Essentially the two states are still located near the 
origin and the surface separately even in the presence of 
the Zeeman field. Also a weak and radial potential V(r) 
can also remove the degeneracy of the two bound states, 
although the two states are still robust. 

In short the exact solutions of a magnetic monopole 
with strength cf>Q in a topological insulator do not fully 
support the Witten effect, a key feature of a axion me- 
dia. This raises a question whether a topological insula- 
tor is really an axion media as predicted by topological 
field theory or not. The finite size effect may produce 
the Witten effect, which suggests that topological insu- 
lators are very close to an axion media. The solutions 
also demonstrate that a magnetic monopole can induce 
spinlcss bound state, which is a rare example of the spin- 
charge separation in three dimensions. 
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